ON MEASUREMENTS OF THE CANONICAL PHASE OBSERVABLE 
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Abstract. Measurements of single- mode phase observables are studied in the spirit of the 
quantum theory of measurement. We determine the minimal measurement models of phase 
observables and consider methods of measuring such observables by using a double homodyne 
detector. We show that, in principle, the canonical phase distribution of the signal state can 
be measured via double homodyne detection by first processing the state using a two-mode 
unitary channel. 



I. Introduction 

In quantum optics, the concept of phase for a single mode electromagnetic field has remained 
a somewhat controversial topic. Alternative descriptions for phase observables have been de- 
veloped, and hundreds of articles (see, e.g., the reference lists of j5UJ E2J GSJ ETJ EH]) and several 
monographs (e.g., |TTj [37]) have been written on the subject since Dirac's famous paper [9j 
published in 1927. A major reason for this variety of phase theories is that in trying to de- 
fine the phase of a quantum oscillator, one can clearly see the restrictions of the conventional 
approach which identifies observables with selfadjoint operators, or equivalently, their spectral 
measures. In fact, it is well known that no spectral measure satisfies all physically relevant 
conditions posed on phase observables (see, e.g., [T%| [55]). and this problem has been over- 
come with the introduction of the more general concept of observables as (normalized) positive 
operator valued measures (POVMs). 

A natural requirement for the description of a phase measurement is covariance with respect 
to phase shifts. In other words, the application of a phase shifter on the field prior to the 
measurement should only shift the phase distribution without changing its shape. Although 
there exist an infinite number of phase shift covariant POVMs, it is generally accepted that the 
canonical phase measurement for the single-mode radiation field is represented by the London 
phase distribution [28j. Hence, the canonical phase measurement is described by the canonical 
phase observable E can , 
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where X C [0, 2ir), {\m)} is the number basis, and \6) = Ylm=o e im6 '|?n) is the (formal) Susskind- 
Glogower phase state [39]. We also recall that the canonical phase measurement arises as the 
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limiting distribution of the Pegg-Barnett formalism [22] ■ In addition, E can has been inde- 
pendently derived by Helstrom [TB] and Holevo [18J in the more general context of quantum 
estimation theory. 

The canonical phase has a number of properties which makes it an optimal choice among 
other phase observables: For instance, E can is pure, i.e. an extreme point of the convex set of all 
POVMs [15]. Any other covariant phase observable F is connected to the canonical phase via a 
quantum channel $ as F(X) = $* (E can (X)) [14] . The canonical phase is also (essentially) the 
only covariant phase observable which generates number shifts [24] . Furthermore, E can and the 
photon number N are noncoexistent , probabilistically and value complementary observables 
P>]. Finally, the canonical phase distribution of coherent states \a), a G C, tends to a Dirac 
delta distribution in the classical limit |a| — > oo [24] . A list of further properties of E can can be 
found, e.g., in p. 51 of [33] . 

The problem of finding a suitable realistic measurement model for the canonical phase ob- 
servable is the last big open problem concerning the quantum description of the phase of an 
electromagnetic field. Of course, E can (or any observable of the field) can be 'measured' in- 
directly or 'sampled' by first measuring some informationally complete observable (e.g. the 
Q-function or the collection of rotated quadratures) and then constructing the canonical phase 
distribution. However, these approaches may hardly be regarded as direct measurements in the 
spirit of the quantum theory of measurement [4j. Some suggestions for direct measurements 
of E can can be found from literature, see e.g. [3J, but they have not led to experiments. It 
seems that from the experimental point of view, the most easily accessible phase observables 
are those arising as the angle margins of certain (translation covariant) phase space observables. 
Such observables are often referred to as (covariant) phase space phase observables, the most 
familiar example being the angle margin of the Husimi Q-function of the field [19j. A natural 
measurement setup for these phase space measurements is then the double homodyne detection 
scheme (also known as eight-port homodyne detection) (see, e.g., [27]), which has also been 
demonstrated experimentally [40]. Even though the phase space phase observables seem like 
a natural choice for measuring the phase of the electromagnetic field, they suffer from certain 
drawbacks when compared with the canonical phase. For instance, the phase space phase ob- 
servables are not pure [6], and the canonical phase E can gives smaller (minimum) variance in 
large amplitude coherent states than phase space phase observables [24] . 

The purpose of this paper is to take a step in the direction of obtaining a realistic measurement 
scheme for the canonical phase. First, we determine those measurement models of covariant 
phase observables which are minimal (in a sense to be defined). In particular, we obtain as 
a special case the minimal measurement models of the canonical phase E can . Since previous 
works indicate that the natural candidates for phase measurements are actually measurements 
of some phase space observables, which would then give the phase distribution after integration 
over the radial part, we then study joint measurements of arbitrary POVMs and the canonical 
phase. After that we focus on the so-called phase shift covariant phase space observables, where 
we also present some physically relevant examples. Finally, we turn to the double homodyne 
detection scheme and show that, in principle, the canonical phase distribution of the signal 
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state can be measured with the setup by adding an extra two-mode unitary coupling in front 
of the apparatus. 

2. Preliminaries 

Throughout this paper we consider a single mode electromagnetic field as our physical system. 
The Hilbert space of the system is therefore % — L 2 (R), and it is spanned by the orthonormal 
basis consisting of the number states \n), n G N = {0,1,2,...}, which in the coordinate 
representation are identified with the Hermite functions. We denote by a*, a and N = a* a 
the creation, annihilation and number operators related to this basis. Moreover, or briefly 
/ is the identity operator of H. We denote by C(H) and T(H) the sets of bounded and 
trace class operators on %. The states of the field are represented by positive trace class 
operators p G Ti^K) with unit trace (i.e., density operators), and the observables are represented 
by normalized positive operator valued measures (POVMs) E : B(Q) — > C(H) on the Borel 
o"-algebra of the topological space Q of possible measurement outcomes^ We say that an 
observable is sharp if it is projection valued, that is, a spectral measure. The measurement 
outcome probabilities are given by the probability measures X i— y tr [pE(X)]. For our purposes, 
the relevant measurement outcome spaces are the torus T, which is identified with [0, 2tt) 
with addition modulo 2tt, the space of nonnegative real numbers M + , and the phase space 
C ~ I 2 ~ 1 + x T. 

While observables represent the statistical description of a measurement, they say nothing 
about the change in the system's state caused by the measurement. The effect of the measure- 
ment on the system is captured in the concept of an instrument [8J. Recall that a completely 
positive instrument is a map X : B(Q) — > £(T("H)), X i-» Xx such that each Xx '■ T(T~l) — > T(H) 
is completely positive, tr[Xx(p)] = 1 for all states p, and the map X i— > Xx(p) is cr-additive 
for all states p. Each instrument determines the dual instrument M. : B(Q) — > C{C{1-C)) by 
tr[pM x {B)\ = tT[X x (p)B] for X G B(Q), p G T{H) and B G C(H), and the associated ob- 
servable E via tr [pE(X)] = tr [Xx(p)] , or equivalently, as E(X) = A4x{I)- in this case we say 
that the instrument (or its dual) is E- compatible. It is well known that, for a given POVM E, 
there are infinitely many E-compatible instruments, that is, an observable can be measured in 
many ways. 

Any completely positive instrument can also be dilated to a quantum measuring process |31j . 
that is, any instrument can be realized as a unitary measurement model of the corresponding 
POVM: For any instrument X there exist a Hilbert space /C associated with the probe system, 
an initial pure state of the probe, a unitary measurement interaction U between the 

system and the probe, and a spectral measure Q of the probe, the pointer observable, which 
define the instrument X uniquely via 

Mp) = tr K [u( P ® |0 (£!)£/*/ ® QPO] 



Note that the domain of definition of a POVM can be a more general a-algebra A consisting of subsets of 
the measurement outcome set f2. 
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where tr^[-] denotes partial trace over the Hilbert space of the probe. Since for each POVM E, 
there exists an infinite number of instruments which are compatible with E, and any instrument 
has infinitely many realizations, it is natural to ask which are the minimal measurement models 
of E. In such a measurement model there would be no unnecessary degrees of freedom in the 
measuring process. This question was solved in |36j . 

We may now turn our attention to the problem of covariant phase observables. As already 
mentioned before, covariance here is with respect to phase shifts which then reflects the re- 
quirement that any phase shift performed on the field prior to the measurement should cause 
a corresponding shift in the distribution of the measurement outcomes. 

Definition 1. An observable E : £>([0,27r)) — > £("H) is a covariant phase observable if 

(1) e i9N E(e)e- i9N = E(& + 6) 

for all G B([0, 2n)) and 9 G [0, 2n), where + denotes addition modulo 2n. 

Condition ([T]) already tells us quite a lot about the structure of the observables. Indeed, the 
following structure theorem for phase observables is well-known, see e.g. |23| . 

Theorem 1. Let E : i3([0,27r)) — > C(H) be a covariant phase observable. Then there exists a 
unique positive semidefinite matrix (c mn ) mjne N with the unit diagonal (i.e. c mm = 1) such that 

oo 

(2) E(0) = ^ c mn — / e i6 ^- n Ue\m)(n\ 

m,n=0 ^ e 

for all G i3([0,27r)). We say that (c mn ) mi „ e N is the phase matrix associated to E. 

Let E be a phase observable with the phase matrix (c mn ) mj „ g N. It is a standard result (i.e. a 
minimal Kolmogorov decomposition of a positive semidefinite matrix) that one can write, for 
all m, n G N, 

where the unit vectors r] m G /C C "H span a Hilbert subspace JC of We say that the ?7 m 's 
are (minimal) structure vectors of E. 

From Theorem 1 of [36] follows that any E-compatible completely positive dual instrument 
is of the form 

oo „ 

M e (B)= V — / ( Vm \T e (B) Vn )e^ m - n Ue\m)(n\ 



Recall that JC and vectors r] m are unique up to a unitary map, i.e. if {r]' m } C JC' is another family of vectors 
satisfying the above requirements then there exists a unitary map U : JC — > KJ such that r( m = U r? m , m € N. 
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for all G £?f [0,27r)) and B G C(H), where any T s : C(Ji) — > C(JC) is a completely positive 
dual channeljj Indeed, by using the (minimal) Kraus decompositions 

r(9) 

k=i 

where A k (9) : JC — > H are bounded operators (such that J2k^k(@)*&k(@) = Ik) an d r{6) G 
N U {00} is the rank of the channel Tq, one sees that 

r{9) 



B^M @ (B)= [ X>' fe (0)*5A' fe (0)f 



' e k=i 

is completely positive; here N k {9) := Afc(#) e~ i£ ' n |?7 n )(n| is a (minimal) pointwise Kraus 
operator of Ai [35J. Moreover, Aie{In) = E(G) so that Ai is E-compatible. We say that Ai is 
rank-1 if r(0) = 1. 

Following Davies [7], we say that an E-combatible instrument Ai is (phase shift) covariant if 
e !W M e (e-' w Be !W )e- jW = M e+e (B) e WN l @ (e-' ld N pe im )e- im = X 0+ ,(p). 

It is easy to check that this holds exactly when 

T g (B) =T(e- i0N Be WN ) 

where T : £("H) — > £(/C) is a completely positive channel. By using a (minimal) Kraus 
decomposition of T, i.e. T(B) = J2l=iK BA k, one sees that T g (B) = YI k =i k k{9)* B k k (6) 
where A k (d) := e i9N A k for (almost) all 9. 

3. Minimal pure measurement models of covariant phase observables 

3.1. General case. We will next consider minimal (pure) measurement models of phase ob- 
servables. In such models, the ancillary or probe Hilbert space is the 'smallest' possible. It is 
quite obvious that any minimal measurement model of E is given by a rank-1 instrument, that 
is, Ai is of the form 

00 

M e (B)= ^ / (Vm\m*BA(6)Vn)e^ m - n) dff\m)(n\ 

m,n=0 

where A(8) : /C — > "H is an isometry for (almost) all 8, see |36| Section 5]. Now the ancillary 
space is (unitarily equivalent with) H! = L 2 ([0,27i~)), consisting of square integrable 'wave 
functions' ip : [0, 27r) — > C, and the pointer observable is the canonical spectral measure 
Q: B([0,27r)) 

(Q(e)VO(0):=Xe(0)VW 



3 Recall that 9 h-> (r]\Tg(B)rf) is measurable for all 77 € CQC), B € and the family {Te} eg [ 2T ) is unique 

(almost everywhere) . 
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where xe is the characteristic function of [36, Remark 3]. Since the vectors e n (8) 
e~ l6n / y/2n, n G Z, form an orthonormal basis of H' one can write 

Q(0)= V (e m |Q(9)e n )|e m )(e n | = V — / n > d9 \e m )(e n \. 



m,n=—oo m,n=— oo 



If the initial pure state £ G H', ||£|| = 1, of the apparatus is given, then one can construct a 
unitary map U : "H ® W — >■ "H (g> "H' by extending 

[^(|n)®£)](0) = (A(% n )«)e n (0). 

Indeed, for all m, n G N, 

tv[U{\n){m\®\i){i\)U*{I u ®Q{Q))} = (U(\m) ® 01(2* ® Q(©))^(k) ® 0) 



([l/(|m>®0]W|[l7(|n>®0]W>de 
^<A(%JA(% n )e^-^ 

(r7mhn>e ifl(m - n) — = tr[|n)(m|E(e)] 

so that, for any state p G T("H), the phase measurement probability distribution i— > tr [pE(0)] 
can be obtained by measuring the pointer Q in the state U(p ® where £ is the initial 

pure state of the apparatus and U describes the measurement interaction between the system 
and probe. 

If, additionally, we assume that the rank-1 instrument (or the minimal measurement) is phase 
shift covariant then one must have A(0) = e iBN A where A : fC -> H is an isometry. In this case, 

oo „ 

M e (B) = M%{B) := V — (Ar] m \e~ im Be im Ar] n ) / e ie{m ~ n) d9\m)(n\. 

— ' 27T Jp, 

m,n=0 J u 



Since c mn = (i] m \i] n ) = (Ar] m \Ai] n ) the vectors Arj m are new structure vectors of E and it is 
enough to study the case A = 1 where lrj := rj for all rj G K. The predual instrument of A4 1 is 



oo „ 



i6N ^i8(m—n) 



/pi 

m,n=0 

and it describes the minimal covariant measurement of the phase observable E with structure 
vectors r\ m . 

Remark 1. We have seen that the minimal ancillary space of any phase measurement is 
(unitarily equivalent with) L 2 ([0,27r)) and the pointer observable can be identified with the 
canonical spectral measure Q. However, the physical interpretation of Q in the case of a single- 
mode optical field is more or less unclear. 

If, instead of the minimal space L 2 ([0,27r)), we consider an ancillary space consisting of two 
modes with the lowering operators a®I% and <g> a (operating on the two-mode Hilbert space 
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% ® H) then Q could be viewed as a spectral measure acting on a subspace of % <S> "H in the 
following way: Choose an isometry J : L 2 ([0, 2tt)) H<S>H. Thus, := JL 2 ([0, 2tt)) C H®H 
is a subspace and the vectors \m)) := Je m , m G Z, form its orthonomal basis. Now 



oo „ 

JQ(0)J* = V — / e' e(m - n) ^|m))((n| 



m,n=— oo 



is a spectral measure on the minimal ancillary space %" . 

For example, one may choose |m)) = \m, 0)) where, for all m G Z, A; G N, 



|m, fc)) :- 



|fc + m) ® | A;), m > 0, 
® |fc — m), m < 0, 



are Ban's relative number states (RNSs) [lj which are the eigenstates of the number difference 
AN = N<g>I-I<g)N and the total number EN = N®I + I®N operators: 

AN\m, k)) = m\m, k)), EiV|m, k)) = (2k + |m|)|m, k)). 

Moreover, they form an orthonormal basis of H <E> H since 

(3) \m) ® \n) = \m — n, min{m, n})), m,nGN. 

Now JQJ* is the restriction of the spectral measure O ang i e : B([0, 2ir)) — > C(H ® H), 

OO „ oo 

(4) anglc (6):= Yl ^Je ie ^deJ2\^,k))((n,k\ 

m,n=—oo k=0 

to H" (see, Remark E]). Note that tr [{p ® |0)(0|)O ang i e (©)] = tr [pE can (0)] where E can is the 
canonical phase. 

3.2. Canonical phase. We will now consider the canonical phase observable E can which is 
determined by the constant phase matrix c mn = 1 or by structure vectors rj m = rj where rj G H 
is any fixed unit vector, that is, 



oo „ 

E can (e)= V — / e^ m - n Ue\m)(n\ 

m,n=0 JIJ 



Now /C = Cf] ~ C and the Kraus operators A(0) : Crj — > H (related to a minimal measurement) 
can be identified with the unit vectors ipg := A(8)r] G "H. Now the corresponding instrument is 

oo „ 

M e {B)= Y 2~ / (MB^)e ie{m - n) de\m){n\ 

m,n=0 ^ e 

and 

[U(\n)®O](0) = ^®e~ Wn . 
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Especially, the phase shift covariant (minimal) instruments are of the form@ 

Ze(p)= E WpWit I j m \r]){r)\e- im j e{m ~ n) de= f e*»<7j| e -*"tr [pE C!m (de)} 

and the related measurement interactions are given by the formula <S> £)](&) = (e l6N rf) ® 

e -i9n ^y^j^ can k e extended to a unitary operator in infinitely many ways. Especially, if one 
chooses rj = |0) and £ = e$ then U(\n) ®eo) = |0) ® e n , n e N, which extends, e.g., to a unitary 
operator [/ : U ® W -»■ ft ® ft', 

— 1 oo oo oo 

tf = E E' n ^ n l® l e fc)( e fcl + EEI^ ra l ® l e «)( e ^l 

A;=— oo n=0 fc=0 n=0 

whose eigenvectors are \n) ® e&, n > 0, < (related to the eigenvalue e l ° = 1) and, for all 
k, n > 0, 

A = 2- 1/2 (|n) ®e fc ± ®e n ), fc ^ n, 
t = |n> ® e n , 

which are related to the eigenvalues +1 = e*° and —1 = e l7T . Hence, U = e lH where the 
generator (Hamiltonian) H of the measurement interaction is 

By using the isometry J of RemarkHJ Je m = |m)) = \m, 0)) = |m) <8> |0) G ft <8>ft when m > 0, 
one sees that 

(In ® J)i)~ k = 2~ 1/2 (|n) <g> Je fc - <g> Je n ) = ^ ® |0) e ft ® ft ® ft. 
where ip~ k = 2~ 1 / 2 (|n) <g> \k) — \k) ® |n)) so that 

ffj := (I n ® J)H(I H ® J)* = n fcX^I ® |0> <0|. 

Since J* J = 1^/, 

£7, := (I n ® J)U(I H ® J)* = (In ® J)e tH (I n ® J)* = e lHj 
and Uj just flips the vectors of the first two modes: 

Uj(ip®(p® |0» = (p®ip® |0), if), (pen. 

^Recall that any E can -compatible instrument is nuclear, i.e. of the form I@(p) = /@ erg tr [pE can (d#)] where 
{°"0}ee[o,27T) C T(H) is a (measurable) family of states [36, Theorem 2]. Now erg can be interpreted as a posterior 
state when the outcome 9 is obtained. It does not depend on the initial state p. In the case of a minimal covariant 
measurement of the canonical phase, ag — e l9N \ri)(ri\e~ 10N for some unit vector 77. The choices 77 = \n), n £ N, 
give the only rotation invariant states <jg = \n)(n\. 
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However, it is still an open question how to measure the pointer observable 

oo „ 

JQ(Q)J* = Y — / e ie{m ~ n) dO\m))Un\. 
2tt J e 

m,n=—oo w 

We will come back to this problem in Remark [3] below. 

4. Joint measurements including the canonical phase 

As we saw in the preceding section, to measure a phase observable E, one needs a pointer 
observable Q a : B(Q) — > £("H a ), Q C R", which satisfies the following requirements: a) the 
Hilbert space "H a of the apparatus must be 'larger' than L 2 [0,27r), i.e. there exists an isometry 
J : L 2 [0,2tt) — > Ha, and b) there exists a pointer function (measurable surjection) / : Q — > 
[0, 2ir) such that 

(5) J*Q a (r 1 (0))J = Q(6). 

For example, if one chooses Q = C ~ R 2 and uses polar coordinates (r, 9) then the choices 
U a = L 2 (C) ~ L 2 (R + ) ® L 2 [0,2tt), f(r,6) :=9, (Q a (Z)*)(z) = Xz(z)*(z), Z G B(C), 
\1/ G "H a , z e C, and (Je n )(r,6) = g(r)e n (9), J °° \g{r)\ 2 dr = 1 satisfy ([5]). Similarly, one can 
choose Q = R 3 (or any R n ) and change the cartesian coordinates to (generalized) spherical 
coordinates. Now one of the coordinates is 6 G [0, 2tt) and we may generalize easily the above 
construction. Hence, an interesting special case of ([5]) is the case of Q = Q' x [0, 2tt) where 
Cl' C R n_1 and f(co,9) = 9 for all 9 G [0, 2ir) and cu G . Now Q a can be viewed as a joint 
measurement of the POVMs ^ Q a (fi' x 6)) and B(Q') 9l4 Q a (X x [0,2vr)). Note that 
if U is the measurement interaction and £ the initial pure state of the apparatus then 

(6) E(0) = V£U*(I n g) Q a (Q' x Q))UV^ 

where : "H — > "H <8> "H a; ^ ^ ® is an isometry. By defining another POVM 

F(X) :=l/ ? T(/,®Q a (Ix [0,2tt))UVz 

we see that E and F are jointly measurable and Q a (resp. U) can be used as a pointer observable 
(resp. measurement interaction) also for F. Next we give the definition for jointly measurable 
observables and study the case of the canonical phase. 

Definition 2. Let E x : B(Vti) — > C{V) and E 2 : B(Q 2 ) — > C(H) be observables. We say that E t 
and E 2 are jointly measurable if there exists an observable E : B(Qi x Q 2 ) ~ > C,(T-L) such that 

(7) Expf) = E(X x 2 ), E 2 (F) = E(fi! x Y), 

for all X G B(fii) and F G B(fl 2 ). 

It is a standard result that for sharp observables joint measurability is equivalent to commu- 
tativity. The same is true also in the case that one (but not both) of the observables is merely 
a POVM. In such a case, E(X x Y) = E 1 (X)E 2 (V) = E 2 (V)E 1 (X). We will see (Prop. □ below) 
that in the case of the canonical phase this rules out any joint measurements where the other 
observable is sharp. This is a consequence of the following lemma. 
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Lemma 1. Let A £ C(U). Then [A, E can (9)] = for all 6 £ B([0,2n)) if and only if A = cl 
for some c £ C. 

Proof. First notice that A can be decomposed as A = A\ + iAi where Aj £ C(7V) are self- 
adjoint; for example, A\ = (A + A*)/2. Therefore, if A commutes with all E can (X), then by 
[A*, E can (X)] = —[A, E can (X)]*, so does A* and it follows that both A,:s also commute with all 
E can (X). Now for any k £ Z define the operator 

B k = — e ik6 dE C3Il {9) = ^ / e^ +m -^^|m)(n| = ^ |m)(m + fc|. 

m,n=0 m=max{0,-/c} 

If commutes with E can , it also commutes with all B k :s. Thus, we have (n|[A,-, -Bfc]|Z) = for 
all n, I £ N and k £ Z. In particular, if k > I then 

(0\[A j ,B k ]\l) = -(k\A j \l) = 

so that Aj must be diagonal. (Recall that Aj is selfadjoint.) But then we also have 

(0\[Aj,B k }\k) = (0\Aj\0)-(k\A 3 \k)=0 

for all k > so that = Cj/ for some Cj £ R Thus, A = cl where c = C\ + ic2 £ C. □ 

Proposition 1. T/ie canonical 'phase observable is not jointly measurable with any nontrivial 
sharp observable. 

Proof. If, say, F : B(Q) — > C(T-L) is a sharp observable which is jointly measurable with the 
canonical phase E can , then F and E can must commute according to |29l Theorem 1.3.1]. By 
Lemma [H (and the fact that F is an observable), for each X £ B(Q) there exists a cx £ [0, 1] 
such that F(X) = Cxi- Thus, F must be trivial. (Moreover, since F is projection valued, one 
must have {cx) 2 = Cx so that X y cx £ {0, 1} is an extremal probability measure which is 
usually a Dirac measure). □ 

The next proposition shows that even if we allow the other observable to be a POVM, we 
still do not gain anything by performing a joint measurement. Indeed, we see that the other 
observable is necessarily a post-processing of E can and therefore it does not give any more 
information than the measurement of merely E can . 

Proposition 2. Let F : B(Q) — > £(%) be an observable. Then F is jointly measurable with 
the canonical phase observable E can if and only if there exists a (weak) Markov kernel m : 
B(Q) x [0, 2tt) [0, 1] such that 

(8) F(X)= m(X,6)dE C!m (d) 

Jo 

for all X £ B(Q). In such a case, the joint observable M : B(Q x [0,27r)) — > is unique 

and is given by 

00 „ 

(9) M(Xx0)= J2 — J_rn(X,6)e ie{m ~ n) d6\m){n\ 



m,n=0 



'e 
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for all X G B(Q) and 6 G B([0, 2tt)). 



Proof. Assume first that F and E can are jointly measurable, with the joint observable M, i.e., 
Mi = F and M2 = E can . Consider the minimal Naimark dilation (/C, Q, V) of E can where 
/C = L 2 ([0,27r)), Q is the canonical spectral measure on /C, and V : H. — > K, is the isometry 
V\n) = e n where e n {9) = ^e~ mi '. Since M(Ix9)<M(0x9) = E can (6) for all X G 
and G i5( [0, 2vr)), there exists an observable R : B(Q) — > C(IC) which commutes with Q and 

(10) M{X x 6) = V*R(X)Q(G)V 

for all X G B(Q) and G -BQ0, 27r)) [T3| Lemma 4.1]. But the canonical spectral measure Q is 
maximal in the sense that any operator which commutes with Q must be a function of it [TU] 
Theorem 1, p. 187]. Indeed, for a fixed X G B(Q) and for any G i3( [0, 27r)) we have 

[R(X) Xe ] (9) = [R(X)Q(6)1] (9) = [Q(6)R(X)1] (9) = [R(X)1] (9) X e(9) 

for almost all 9 G [0,27r) where 1 denotes the constant function 9 1— > 1. Since the linear 
combinations of characteristic functions are dense in /C, we conclude that R(X) corresponds to 
multiplication by the function R(X)1. By denoting m(X, 9) = (R(X)1)(#) we obtain a (weak) 
Markov kernel m : B(Q) x [0, 2n) — > [0, 1], and it follows from Eq. f lTII]) that 

F (e) = F*R(X)y= / m(X,^)rfE can (0). 

Suppose now that there exists a (weak) Markov kernel such that Eq. ([8]) holds. Since for each 
X G i3(f2) the map 9 i-> m(X, ^) is measurable and m(X, 0) < 1 for almost all we may define 
a bounded operator R(X) G £(/C) by (R(X)<^)(^) = m(X,9)<p(9). The map X h-> R(X) is then 
a POVM which commutes with Q. Therefore the map (X, G) 1— > V*R(X)Q(<d)V extends to an 
observable M : B(Q x [0,27r)) — > £("H) whose margins are F and E can . In other words, F and 
E can are jointly measurable. 

In both of the above instances the joint observable M satisfies Eq. (Q, and the uniqueness fol- 
lows from [T3| Theorem 4.1(a)] since E can is an extreme point of the convex set of all observables 
on S([0,2tt)) [15]. □ 

Note that, by combining the preceding two propositions, we see that F is projection valued 
if and only if the (weak) Markov kernel m(X,9) = Cx = {cx) 2 (almost all 9), i.e. it does not 
depend on 9 and cx G {0, 1}. 

5. Phase shift covariant phase space observables 

Even though we have already seen that no additional information can be obtained by mea- 
suring some other observable together with the canonical phase, there is a strong practical 
motivation for studying such joint measurement. For instance, the realization of the measure- 
ment of the canonical spectral measure Q : i3([0,27r)) — > L 2 ([0,27r)) may be problematic, but 
things get much simpler if we instead consider Q as the angle margin of the canonical spectral 
measure M : £>(C) — > L 2 (C) ~ L 2 (IR + ) <g> L 2 ([0,27r)), i.e., we interpret Q as the observable 
Q(^0 011 the larger Hilbert space. In this case L 2 (C) may be simply realized as the 
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Hilbert space of a two-mode field, whereas M, and therefore Q, can be measured by performing 
homodyne detection on the two modes. 

For any observable P : B(C) — > C(H) we define the angle and radial margins P an gie : 
B{[0, 2tt)) C{U) and P rad : B(R+) -> C{U) via 

(11) Pan g ic(e) = P(R + x 9), P rad (X) = P(X x [0, 2tt)). 

Any phase observable E can obviously be obtained as the angle margin of a phase space ob- 
servable in a trivial manner. Namely, for any probability measure fi : £>(M+) — > [0,1], the 
observable P : B(C) ->■ C(H) for which P(X x 6) = /i(X)E(6) is a suitable choice. This 
observable has the additional property that it also transforms covariantly under the action of 
the phase shifter, i.e., 

e i9N P(X x Q)e~ im = /j(X)e ieN E(Q)e~ ieN = /i(X)E(0+0) = P(X x (6+0)). 

We will now focus our attention on phase space observables which have this symmetry property. 

Definition 3. A phase shift covariant phase space observable is a POVM P : B(C) — > £(TL) 
such that 

e i9N P(X x Q)e- ieN = P(Xx (6+0)) 
for allXe B(R+), 6 E B([0,2tt)), and 6 E [0,2tt). 

The following structure theorem for phase shift covariant phase space observables can then 
be proved. 

Proposition 3. An observable P : B(C) — > £(91) is a phase shift covariant phase space 
observable if and only if there exists a probability measure fi : B(R + ) — > [0, 1] and (a weakly 
^-measurable field of) vectors t] m (x) G % satisfying the condition 



\rj m {x)\\ 2 dfi(x) 



such that 



(12) P (Xx0)= V — / (e- im V(^)|e- m V(^))^M^)l^)(^l 

z — ' Z7T IxvPi 



m,n=0 

for allXeB(R + ) and 6 G B([Q, 2tt)). 



Proof. Consider a fixed set X G B(R+). Since e ieN P(X x e)e- ieAr = P(X x (6+0)) for all 
G B([0, 27r)), following the proof of the structure theorem for covariant phase observables 
|23| . there exists a positive semidefinite matrix (c mn (X)) mngN such that 



oo „ 

P(X x 0) = c m n(X)— / e« m ~ n ^d9\m)(n\ 



m,n=0 
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Now, for all m,n GN the map X \— > c mn (X) is a complex measure which is clearly absolutely 
continuous with respect to the probability measure 

oo 

X ^ fi(X) = Y^ =0 ^n(n\P{X x [0,27r))|n> = ^A n c nn (X) 

n=0 

where A n > for all n and X n — 1. It follows from [20J that there exist vectors r] m (x) G % 
such that c mn (X) = j x (r] m (x)\r] n (x))dfj,(x). Hence, Eq. (fT2l) holds. The converse claim is 
clearly true. □ 

Using Prop. [3] we can now determine the angle and radial margins of any phase shift covariant 
phase space observable. Indeed, by setting X = M + in Eq. (fl2"|) we obtain 

oo 

(13) Pan g ie(e) = £ c mn — / e"^ n ^ d9\m) (n\ 

m,n=0 ^ 

where the phase matrix elements are given by 



(r) m (x)\r) n (x))dn(x). 
Similarly, the radial margin is seen to be 

oo 

(14) P Iad (X) = Y J Cmm{X)\m)(m\. 

m=0 

The map (X, m) (->■ c mm (X) is a Markov kernel on i3(R + ) x N, and the observable P ra d can 
therefore be viewed as a smeared number observable. More generally, we say that an observable 
F : B(Q) C(jH) is a smeared number observable if there exists a Markov kernel m : B(Q) x 
N -> [0, 1] such that 

oo 

F(X) = J2m(X,n)\n)(n\ 

n=0 

for all X G 

Proposition 4. Let F : <B(R+) — > 6e a smeared number observable and let E : £>(IR + ) — )■ 

£("H) &e a covariant phase observable. Then E anc? F /lave a joint observable if and only if they 
have a joint observable which is a phase shift covariant phase space observable. 

Proof. Assume that E and F have a joint observable M : B(C) — > £("H) and define the biob- 
servable P' : B(R+) x B([0, 2vr)) C{U) via 

(15) P'(X,6) = — / e~ ieN M(X x (0+6))e ieN d6 

where the integral is understood in the weak*-sense. By [26| Theorem 6.1.5], P' extends to a 
unique observable P : B(C) — > £■{%) which is clearly phase shift covariant. Setting O = [0,27r) 
(resp. X = 1R + ) in Eq. (fT5l) and using the phase shift invariance of F (resp. phase shift covariance 
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of E) we then see that P ra d = F (resp. P an gie — E). Hence, E and F have a joint observable which 
is a phase shift covariant phase space observable. The converse statement is trivial. □ 

We immediately see that the canonical phase observable is obtained as the angle margin of 
a phase shift covariant phase space observable if and only if 

/■oo 

(Vm\Vn) ■= / (i] m {x)\r] n {x))dfi(x) = 1 
Jo 

for all m,n G N. In particular, we must have (//ol^n) = 1 for all n G N and since rj n is a unit 
vector, this implies that r\ n = rj for all n G N. But in such a case the radial margin is the 
trivial observable P ra d(^) = Cqq{X)I. A similar situation happens if we insist that the radial 
margin is the sharp number observable. In this case we have c mm (X) = 5 m (X) (the Dirac 
measure concentrated on m) so that, in particular, J, , \\r] m (x)\\ 2 dfi(x) = 1 and we must have 
Vm(x) = for /i-almost all x G M + \ {m}. It follows that c mn = {i] m \r] n ) = 5 mn and the angle 
margin is therefore the trivial one: P an gie(@) = ^ le d6 1. 

Example 1. As a first example, we consider phase space observables which, in addition to 
phase shift covariance, are covariant with respect to phase space translations represented by 
the displacement operators D(z) = e za *~ za , z G C. Generally, any phase space POVM G which 
satisfies the covariance condition D(z)G(Z)D(z)* = G(Z + z), z G C, Z G B(C), is generated by 
a unique positive trace one operator T, giving the corresponding observable the explicit form 

(16) G T (Z) = - [ D(z)TD(z)*d 2 z. 

Jz 

However, not all of these observables are covariant with respect to phase shifts. Indeed, this is 
the case if and only if the generating operator is a mixture of number states, T = YlT=o ^k\k)(k\, 
see [23J for the proof of the 'only if -part and the calculation below. In order to connect 
these observables to the general structure discussed above, let us consider G' fe \ the observable 
generated by a single number state First note that we may write Eq. (fTBT) as 

00 i r 

(17) G |fc) (^ x 6) = V — / (m\D(re ie )\k)(n\D(re ie )\k)dr 2 d6\m)(n\ 

where the matrix elements of the displacement operators are given by 

(18) (m\D(re W )\k) = (-l)™x{0,fc-m> f^^m e ^(m-k) \m-k\jjm-k\ ( r 2) e -4 
V / \ I \ / 1 / v / y (max{m,A;})! mm{m,fc}\ / 

and is the associated Laguerre polynomial (which generalizes the Laguerre polynomials 
L n = L°). We can then define the probability measure /i : B(R+) [0, 1] via dfi(r) = e~ r dr 2 , 
as well as the vectors 

n k I \ _ / i Nmax{0,fe-m} / (min{m,fc})~ |m-fc| r \m-k\ , 2\ p nj 

Vm[ r ) - { l ) V (max{m,fc})! r ^min{m,fc}V > ^k G «. 
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where (p^ G % is an arbitrary fixed unit vector. The normalization condition J ||?7m( r )ll 2 dfi(r) = 
1 is then satisfied, and it is merely a simple observation that G' fc ^ takes the form of Eq. (TT2]) . 
Note that the angle margin G^ gle is never the canonical phase observable [23| |H] . 

Example 2. A second example can be obtained by modifying the above considerations. In- 
deed, if we replace the vectors 77^ (r) by £m( r ) — -^min{m,fc}('" 2 ) ( ^fc while keeping the measure ^ 
unchanged, we get an observable F' fc ^ : B(C) — > C(H) with the explicit form 

00 

(19) Fl fc >(X x 0) = J~ / e ie ( m - n )L min{m , fc} (r 2 )L min{ „, fc} (r 2 )e- r2 rfr 2 ^|m)(n|. 

Unlike the observable G' fc \ this observable is not translation covariant, though it is clearly 
covariant with respect to phase shifts. The significant features of these observables are in their 
margins. For the angle margin phase observable Fj^ gle , we get the phase matrix elements 



00 

imin/m. k\ (t )Lmir<{n k\ (T ) 6 df = 5 n 



■'min{m,k}y J-^min{n,fc} y "min{m,fc},min{n,fc} • 

In particular, the angle margin of F' ^ is the canonical phase. On the contrary, if we increase 
the value of k we find that the observable becomes in some sense more trivial. Indeed, since 
c-mn — $mn for m,n < k, the phase matrix always contains a (k + 1) x (k + 1) identity matrix 
in the upper left corner. As for the radial margin, we get 

k oo 

^(x) = ^2c mm (x)\m)(m\ + Ckk (x) Yl ' m >H 

m=0 m=k+l 

which shows that by increasing the value of k the observable becomes, vaguely speaking, more 
and more nontrivial. We also immediately recognize the obvious fact that for k = correspond- 
ing to the angle margin being the canonical phase, the radial margin is trivial. 

Example 3. As a final example, we consider the case where the probability measure is a Dirac 
measure fi = 5 ro concentrated at some point r > 0. The corresponding phase space observable 
P r ° : B(C) — > £(H) can then be written as 

oo i r 

(20) P r °(X x 9) = (VmK) S ro (X)— / e ie ^ d9 \m)(n\. 

m,n=0 ^ e 

We now immediately see that the angle margin P^ie * s the phase observable with the phase ma- 
trix elements c mn = (i] m \r] n ), and the radial margin is the trivial (sharp) observable P^jPO = 
S ro (X)I. 

6. Double homodyne detection scheme 

We will now turn our attention to the double homodyne detector (also called the eight-port 
homodyne detector), a well established method for measuring translation covariant phase space 
observables related to a single mode electromagnetic field [27J[2lJ[25]. This scheme makes use of 
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the fact that the usual single homodyne detector provides a measurement of an arbitrary field 
quadrature. Indeed, when the signal field under investigation is coupled to a strong auxiliary 
field in a coherent state \z) using a 50 : 50 beam splitter, and the scaled photon number 
difference A ( I <g> iV — N £g> /) of the two output modes is measured, then for a sufficiently large 
\z\ the measured observable is approximately the quadrature observable Qg, where 9 = axgz 
and Qe(X) = e %dN Q(X)e~ t9N , Q being the canonical spectral measure on the real line (i.e., the 
position observable) (for more details, see, e.g., |22|). 

In double homodyne detection, the signal field is first coupled to a parameter field in some 
state cr via a 50 : 50 beam splitter, after which a phase shift of — | is performed on one 
of the output modes (see Fig. [T|). Balanced homodyne detection is then performed on each 
output mode, so that by choosing the phase of the auxiliary coherent field to be zero in both 
measurements, this corresponds to measuring the canonical spectral measure on M? for the two- 
mode field. With this setup, the measured observable is the phase space observable G CT where 
the generating operator is connected to the state of the parameter field via the conjugation 
map {C<p)(x) = <f(x) as a' = CaC [21]. In particular, if we want the observable to be phase 
shift covariant and thus give rise to a phase observable as the angle margin, we must use a 
parameter field which is diagonal in the number state representation, a = Y2k=o ^k\k)(k\. The 
simplest case is obviously obtained by using the vacuum a = |0)(0|. 




Figure 1. Schematic of a double homodyne detector. The signal field is mixed 
with a parameter field by means of a 50 : 50 beam splitter, after which a phase 
shift of — 7r/2 is performed on one of the modes. Balanced homodyne detection 
is then performed on both modes. 

The problem with using this measurement setup is of course caused by the fact that the 
canonical phase cannot be obtained as the margin of any translation covariant phase space 
observable [23, 24J. Thus, a modification of the setup is needed. We will next show that a 
suitable modification is obtained by adding a unitary coupling between the signal and parameter 
fields prior to the beam splitter. In other words, while in the usual double homodyne detection 
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the signal and parameter fields are uncorrelated before entering the beam splitter, with this 
modification they will enter the beam splitter in an entangled state. With this method, it is in 
principle possible to measure also the observables F' fc ^ encountered in Example [21 

To this end, first notice that in the usual double homodyne detector, the overall unitary 
coupling consisting of the beam splitter and the phase shifter is given by U : "H; n — > Hout, 

(21) (U\m) <g> \n)) (r, 6) = -= (n\D{re ie )*\m), 



where TL m ~ % (g) H consists of the input signal and parametric field modes and % ou t — 'H®'H 
is the output space. The observable measured with the two homodyne detectors is then (or 
rather may be chosen to be) the canonical spectral measure Q : B(C) — > £("H ut)- If P and a 
are the states of the signal and parameter fields, respectively, then a direct computation shows 
that 

tr [U(p <g> a)U*Q(Z)\ = tr \pG a ' {Z) 

for all Z G B(C). Hence, by considering only the angle margin Q ang ie as the pointer observable, 
one can measure the angle margin G CT '. In particular, by preparing the parameter field in a state 
which is diagonal in the number basis, one can measure any phase space phase observable. 

Now it is easily seen from the discussion in Example [2] that the total coupling needed for the 
measurement of is 

(22) (V\m) ® \n)){r,6) = -^e' 9 ^ L min{mM (r 2 ) e -^ . 

\/7T 



We want to exploit the fact that we already have at our disposal the unitary coupling (|2T|) 
of the usual double homodyne detector. Therefore, we will look for a unitary operator W : 
Hin — > Hin such that V = UW, which would then amount to adding an extra component to 
the measurement setup prior to the beam splitter, see Fig. [2J 
The action of the operator W is now given by 



W(\m) ® \n}) = \ U \ k ) ® 10 V \ m ) ® \ n )) \ k ) ® 10- 

k,l=0 

Let us denote aM,mn = (U\k) <8> |Z) | V"|m) (g) |n)) so that 

a kl , mn = (-l) maxW - fe} yiiI^ I e-(-^--)^-lL^ fci/}( ^)L mm{m , n}( ^)e- 2 rdrdB 

. 1)m ax{0,n- m } / min{fc,fc+n-m}! f°° x \\n-m\ L \n-rn\ ( )L ( x ) e - dx . 



Jl,k+n—m\ 







The above integral can be further calculated, but it seems that a simple closed expression does 
not exist. However, note that this already tells us that H^m) £g> |n)) is an eigenvector of the 
photon number difference operator AiV. Indeed, ANW(\m) £g> |n)) = (m — n)W(\m) ® \n)) 
and Wdm) ® |n)) can be written as a series with respect to the RNSs |m -n,s)),s£ N. 



5 One could of course consider equally well a unitary operator W' : H ut "Hout such that V = W'U, but 
this leads to a similar treatment. 
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Figure 2. A modified double homodyne detector where an additional unitary 
coupling W is performed prior to the beam splitter. By choosing the parameter 
field to be in the vacuum state, on can obtain the canonical phase observable as 
the angle margin of the measured phase space observable. 

Since the canonical phase is obtained as the special case n = of the above consideration, 
we will next focus on that. Now the above integrals are easily calculated and we obtain the 
expressions W(\0) <g> |0)) = |0) <g> |0) and, for m > 0, 

00 TYi* I m I 1 'l 

(23) = ri ± -^au, *))=^h± r<* + f > ( m : k ) " V *» 

where \k + m) <8> \k) is the RNS \m,k)). From V = UW one sees that the canonical phase 
distribution i-> tr[pE can (B)] of any state p can be measured directly by first producing an 
entangled two-mode input state S(p) = W(p <8> |0)(0|)W / * and then measuring S(p) via double 
homodyne detection. 

Remark 2. Note that W(|m) ® |0)) is a two mode nonlinear coherent state |41| which are 
generally of the form \a, f,m) = YlT=oC'k\k + m) (8) \k) where the constants Ck G C are such 
that J^ fc |Cfc| 2 = 1. Moreover, they satisfy the equations 

f(N ® 1,1 ® N)(a® a)\a, f,m) = a\a,f,m), a E C, 

where /(iV <8> I, I <8> iV) is some fixed function of the single- mode number operators. As shown 
in 0U, 




m! 



+ to)! 



1 

II /( s _ i _|_ m> s _ i) 
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so that W(|m)<8)|0)) = \a,f,m) where a — 1, C = T(f + 1)/Vm!, and /(ni, n 2 ) = 2/(m+n 2 ), 
rii, n 2 6 N, fij / 0. 

Remark 3. Note that there does not exist bounded (or unitary) operators D{re l °) such that 



D{re ld )*\m) = e %e<yn m ^ L m i n i m>n \(r 2 )e a for all r and 6*. However, we can write 



(V\m) <g> \n))(r, 9) = -^=((9,r 2 \ \m) <g) \n) ) 

\/7T 



where |0, a;)) := X]^L-oo Sfclo e m6, L fc (x)e _:E//2 |n, k)), 8 G [0, 27r), x > 0, is (essentially) a gener- 
alized vector introduced by Ban [21 Eq. (19b)], see Eq. ([3]). Now the angle margin observable 
of the (two mode) spectral measure 

: B(C) -> C{U®U), Z ^ 0(Z) := V*Q(Z)V = - [ \9,r 2 )) ((9,r 2 \rdrd9 

Jz 

is O ang ie defined in Eq. (jlj), i.e. 0(R + x 6) = ang i e (9). The radial margin is 

OO OO „ 

O rad (X):=O(Xx[0,27r))= L fc (r 2 )L i (r 2 )e- r2 rfr 2 |n,A;))((n,/|. 

n=-oo k,l=0^ X 

Immediately one sees that tr[(p<g> |Jfe)(ife|)0(Z)] = tr [pFl fe >(Z)] . 

Example 4. The action of V is easy to calculate for some physically relevant states. For exam- 
ple, if the first (or the second) mode is in the coherent state \a) = e~' Q ' I 2 Ylm=o o; m / y/m\\m) 
and the second mode is in the vacuum state |0) then 

e -H 2 /2 °^ a m r2 

(V\a) ® |0»(r, 9) = —=- £ " e -*» e - T = -^<0|a), 
that is, we obtain the London distribution 9 i— > (0|a) of the coherent state. Let then 

ip a := C(a) > — rim) ® m), a G C 



m! 

m=0 



be a pair-coherent state where C(a) := Jo(2i|a|) x / 2 is a normalization constant and 

m=0 v y 

is the zeroth Bessel function of the first kind. Now 

(V^ a )(r,9) = ^ E ^ m (r 2 )e- 2 / 2 = ^Jo(2rv^Ke-^ 2 
V7T ' m! V 71 " 

v m=0 
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whereas, for a 'two-mode phase coherent state' ip^ := (1 — |a| 2 ) 1//2 ^^ =0 a m \m) (g) \m + q), 
q G N, a G C, \a\ < 1, one gets the Gaussian function 

~ n I 1 2\ 1/2 00 /1 UJ2U/2 

(y#)(r,e) = iL_^i_ e ^y^ a m L m {r 2 )e- r2 ' 2 = {L ~ J , e ^-[<*/(i-«)+i/2]r 2 
V 71 " n VTr(l-a) 

v m=0 ' 

(see Eqs. 8.975(3) and 8.975(1) of [12]). Finally, we note that, by using V, it is possible to 
produce the (weighted) monomial distributions r 2n e~ r I 2 in the radial direction: Indeed, the 
(nonnormalized) state ip n := J2m=o(~ l) m (m) \ m ) ® \ m ) gi yes 

<y<p n )(r,0) = -|= it(-ir( n )L m (r 2 )e- r2/2 = -^r 2 ^ 2 . 

v m=0 - v 



7. Conclusions 

We have studied measurements of covariant phase observables, with the aim of obtaining a 
realistic measurement model for the canonical phase. To this end, we have first determined the 
minimal measurement models of phase observables. However, since practical issues suggest that 
a realistic phase measurement is most easily obtained by measuring a phase space observable, 
we have considered such joint measurements. As a special case, we have focused on the so- 
called phase shift covariant phase space observables. In particular, we have shown that the 
canonical phase may be obtained as the angle margin of such an observable. Finally, we have 
considered the quantum optical double homodyne detection scheme, and its modification as a 
means of measuring phase space observables. We have constructed a unitary coupling which, 
when placed in front of the setup, would allow one to measure the canonical phase observable. 
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